Three-terminal thermoelectric transport through a molecule 
placed on an Aharonov-Bohm ring 
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The thermoelectric transport through a ring threaded by an Aharonov-Bohm flux, with a molec- 
ular bridge on one ol its arms, is analyzed. The charge carriers also interact with the vibrational 
excitations of that molecule. This nano-system is connected to three terminals: two are electronic 
reservoirs, which supply the charge carriers, and the third is the phonon bath which thermalizes 
the molecular vibrations. Expressions for the transport coefficients, relating all charge and heat 
currents to the temperature and chemical potential differences between the terminals, are derived 
to second order in the electron-vibration coupling. At linear response, all these coefficients obey 
the full Onsager-Casimir relations. When the phonon bath is held at a temperature different from 
those of the electronic reservoirs, a heat current exchanged between the molecular vibrations and 
the charge carriers can be converted into electric and/or heat electronic currents. The related trans- 
port coefficients, which exist only due to the electron-vibration coupling, change sign under the 
interchange between the electronic terminals and the sign change of the magnetic flux. It is also 
demonstrated that the Aharonov-Bohm flux can enhance this type of conversion. 

PACS numbers: 85.65.-|-h,73.63.Kv,65.80.-g 

Keywords: thermoelectric transport, charge and heat transport, electron-vibration interaction, molecular 
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I. INTRODUCTION 

Thermoelectric effects in bulk conductors usually ne- 
cessitate breaking of particle-hole symmetry. In meso- 
scopic structures this asymmetry may be fairly high and 
can be also controlled experimentally, leading to the pos- 
sibility of realizing a relatively large thermopower coef- 
ficient. As a result there is currently much interest in 
investigations of thermoelectric phenomena in nanoscale 
devices at low temperatures. Experimental studies have 
been carried out on point-contactSji^ quantum dots?^^— 
nanotubes,— ^— silicon nanowires,— and more. Very re- 
cently, the paramount importance of another symmetry 
breaking has been pointed out. It has been proposed 
that the thermal efficiency (defined at steady state as 
the ratio of the output power to the heat current) can 
be significantly enhanced once time-reversal symmetry is 
broken and the thermopower coefficient becomes asym- 
metrical (as a function of, e.g., a magnetic field) 

Early theoretical studies of thermoelectric transport 
coefficients of microstructures were based on the Lan- 
dauer approachii"— and were mainly focused on charge 
and heat currents between two electronic terminals, with- 
out coupling to phonons. Later on, effects of electron- 
electron processes and electronic correlations (increas- 
ingly important at lower temperatures), as well as that of 
an applied magnetic field, on the thermopower produced 
in largei^ and single-levelii quantum dots, and also in 
quantum wires^^ were considered. The signature of at- 
tractive electronic interactions on the thermopower was 
considered in Ref. and the dependence of the ther- 
moelectric response on the length of the atomic chain 
connecting the leads has been recently computed within 



a density-functional theory<^ 

The coupling of the charge carriers to vibrational 
modes of the molecule should play a significant role in 
thermoelectric transport through molecular bridges, even 
more so in the nonlinear regime Indeed, a density- 
functional computation of the nonlinear differential con- 
ductance of gold wires attributed changes in the TV char- 
acteristics to phonon heating ;^^:^^ and the thermopower 
coefficient was proposed 0jS 0, tool to monitor the ex- 
citation spectrum of a molecule forming the junction 
between two leadsi ^^'^^ It was suggested that the See- 
beck effect in such bridges can be used for converting 
heat into electric energy;^ and to determine the loca- 
tion of the Fermi level of the transport electrons relative 
to the molecular levels, and also the sign of the dominant 
charge carriers, either for a molecular conductor )22r— or 
for an atomic chaini^i^ This was confirmed experimen- 
tally: the Seebeck coefficient, as measured by STM on the 
benzenedithiol family sandwiched between two gold elec- 
trodes, showed that the charge carriers are holes passing 
through the HOMO, whose location with respect to the 
metal Fermi level was determined from the magnitude of 
the coefficient 1^ 

Theoretically, when the coupling to the vibrational 
modes is ignored, the transport coefficients have the same 
functional form as in bulk conductors, with the energy- 
dependent transmission coefficient and its derivative re- 
placing the conductivity^Si^ Although the corrections to 
the thermoelectric transport due to the coupling to the 
vibrational modes is often small, their study is of interest 
because of fundamental questions related to the symme- 
tries of the conventional transport coefficients, and since 
they give rise to additional coefficients connecting the 



heat transport in-between the electrons and the vibra- 
tional modes. In a recent article^ (referred to below as 
I) wc have analyzed the thermoelectric phenomena in a 
molecular bridge, and studied effects induced by the cou- 
pling of the charge carriers with molecular vibrational 
modes. In particular we have considered the situation 
in which the molecule is strongly coupled to a heat bath 
of its own and thus may be kept at a temperature dif- 
ferent from those of the source and sink of the charge 
carriers, making the junction a mixed thermal/electronic 
three-terminal one. Namely, we have assumed that the 
relaxation time due to the coupling of the molecule to 
its own heat bath, ry, is short on the scale of the cou- 
pling of the molecule to the charge carriers. The latter 
is determined by the electron-vibration coupling, 7 (and 
the conductance electrons density of states) which is the 
small parameter of our theory. Hence, h/r^ may still be 
very small on all other physical scales, such as ftwg, where 
Wg is the frequency of the vibrations, or the molecular 
(electronic) level width. The phonon bath may be real- 
ized by an electronically insulating substrate or a piece 
of such material touching the junction, each held at a 
temperature Ty . Similar experimental three-terminal se- 
tup s, consisting of quantum dots, were discussed in Refs. 
[35L [36I and [ll- The effect of thermal probes on the elec- 
tronic heat conduction have been considered theoretically 
for tunnel normal-superconducting junctions^ and for 
mesoscopic conductorsi^ 

In this mixed thermal/electronic junction, heat sup- 
plied by the phonon baths [to which the electrons and the 
molecule are (separately) coupled] can be exchanged be- 
tween the transport electrons and the vibrational modes. 
Even more intriguing, heat carried from the phonon bath 
to the vibrations can be converted into a charge current 
(or an electronic heat current) flowing between the elec- 
tronic reservoirs. We have found in I that in order for this 
conversion to take place, one needs to break the spatial 
symmetry of a molecular bridge, in addition to the broken 
electron-hole symmetry. Here we examine the outcome 
of breaking another symmetry, that of time-reversal. To 
this end, we exploit the Aharonov-Bohm effect, placing 
the molecular bridge on one arm of a ring threaded by a 
magnetic flux, see Fig. [T] As we show, for a non-zero 
flux the above conversion can also appear in a spatially 
symmetric bridge, provided that the effective couplings 
to the leads are energy dependent. The specific quantities 
we study are the heat and the charge currents in a small 
mesoscopic (or nanometric) system, depicted schemati- 
cally in Fig. [T] a molecule (having one electronic level) 
is placed on an Aharonov-Bohm ring, which is attached 
to two electronic reservoirs (held in general at different 
temperatures, j^, and at different chemical potentials, 
/iL p{ ) . The vibrational modes of the molecule, with which 
the charge carriers may exchange energy, are thermalized 
by a third terminal, kept at the temperature Ty. 

Since the calculations are rather technical, we organize 
the paper as follows. After a short analysis of the entropy 
production in our system (Sec. |lT]), wc summarize the 
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FIG. 1: A mixed, thermal-electronic, three-terminal system, 
modeled by a localized level placed on an Aharonov-Bohm 
ring which is attached to two electronic reservoirs, having dif- 
ferent chemical potentials and temperatures, ajl.r and Tl,r, 
respectively. An electron residing on the level interacts with 
its vibrational modes. The population of these modes can be 
determined by a coupling to a phonon source kept at tempera- 
ture Tv- The ring is threaded by a magnetic flux measured 
in units of the flux quantum. 

main results and in particular discuss the magnetic-flux 
dependence of the transport coefficients. We then pro- 
ceed in Sec. [Hi] to specify our model, and to derive the 
explicit expressions for the currents. Section ITVl details 
the various transport coefficients of the three-terminal 
junction in the linear-response regime. In particular, we 
examine the various symmetries of the coefficients, and 
verify that they obey the Onsager relations. Section |V] 
is devoted to a discussion of our results; we consider the 
transport coefficients in certain simple cases, in partic- 
ular their dependence on the spatial and electron-hole 
symmetries of the junction. 

II. ENTROPY PRODUCTION AND 
THERMOELECTRIC TRANSPORT 
COEFFICIENTS 

The consideration of the entropy production in the 
linear-response regime, carried out in I, is quite illu- 
minating. Using the thermodynamic identity TdS ~ 
dE — ^dN ^ one finds that the entropy production at 
the left (right) electronic reservoir is 

^L(R) = (^L(R) - /^L(R)(^L(R)>) ■ (1) 

Here, —Ei^^yV) is the energy current emerging from the 
left (right) reservoir, while — (iVL(R)) is the correspond- 
ing particle current. Adding to Eqs. ([!]) the entropy pro- 
duction of the vibrational modes, S'v = Ey/Ty, where 
—Ey is the energy current from the phonon bath into the 
molecule, yields the total dissipation of the system, 

+ jr{Ei^- /iL(^L)) + ^ {e^. - t^RiNn)) ■ (2) 
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We demonstrate below that our three-terminal junction 
conserves charge, i.e., (N^ + N^) = 0, and also energy, 
i.e. + + Ey = 0. In the linear-response regime all 
three temperatures (see Fig. [1]) are only slightly different 
from each other, 



^L(R) ~ =^ ~ ' 

Tv = T + ATv , 



(3) 



and the chemical potentials differ by a small amount. 



Ml(r) = M ± ^ • (4) 



One then finds 



Sy + S^ + S^ = ^i-Ey) + ^I+^I<,, (5) 

where / is the net charge current flowing from the left 
electronic reservoir to the right one. 



(6) 



while Iq is the net heat current carried by the electrons. 



Iq = Ie- (M/e)/ , with I^^--i^E^- E^) . (7) 

We keep the notation —Ey for the the energy/heat cur- 
rent from the phonon bath. Thus, the entropy production 
of our three-terminal system is a simple example of the 
general expressions for linear transport, consistent with 
the Onsager theoryJi 

The calculations in the next two sections yield the re- 
lations among the driving fields and the currents. In the 
linear response limit one has 



(8) 



where the matrix of the transport coefficients, A4, is 



/ 




Afi/e 




= M 


AT/T 






. ATV/T _ 



M = 



G($) 
K(-$) 
XV(-$) 



K($) 
XV(-$) 



XV($) 
XV($) 
CV($) 



(9) 



This matrix satisfies the Onsager-Casimir relations. We 
find below that all three diagonal entries are even func- 
tions of the flux $. The off-diagonal entries of consist 
each of a term even in the flux, and another one, odd in 
it, obeying altogether Mtj{^) = Mji{-^). 

The flux dependence of the transport coefficients is 
very interesting. We find (see Sees. |lll] and |IV] for 
details) that there are three types of flux dependencies 
hiding in those six coefficients. First, there is the one 
caused by interference. The interference processes mod- 
ify the self energy of the Green functions pertaining to 
the Aharonov-Bohm ring, in particular the broadening 



of the electronic resonance level due to the coupling with 
the leads. The interference leads to terms involving cos 3>. 
Secondly, there is the flux dependence which appears in 
the form of cos2<I> (or alternatively, sin^ This reflects 
the contributions of time-reversed paths. These two de- 
pendencies are even in the flux. They yield the full flux 
dependence of the diagonal entries of the matrix A4 , Eq. 

and the even (in the flux) parts of the off-diagonal el- 
ements. Finally, there is the odd dependence in the flux, 
that appears as sin4>. This dependence characterizes the 
odd parts of the off-diagonal entries of Ai . These terms 
necessitate the coupling of the electrons to the vibrational 
modes. 



III. THE CURRENTS 

In this section we present our model system, give de- 
tails of the currents' calculation, and verify their balance. 



A. The model 

In our analysis, the molecular bridge is replaced by 
a single localized electronic level, representing the low- 
est available orbital of the molecule; when a transport 
electron resides on the level, it interacts (linearly) with 
Einstein vibrations. Our analysis does not include elec- 
tronic interactions, but focuses on the electron-phonon 
ones. Thus, the Hamiltonian of the molecular bridge, 
which includes the coupling with the vibrations, reads 



- eoCoCo + ^oib^b +-)+ 7(6 + b^ylc^ , (10) 



where eg is the energy of the localized level, Wg is the 
frequency of the harmonic oscillator representing the vi- 
brations, and 7 is its coupling to the transport electrons. 
The Hamiltonian describing the tunneling between the 
molecule and the leads is 

^coup = E(^fe4co + He) + ^(ypctcg + He) (11) 



[using k{p) for the left (right) lead]. The leads' Hamilto- 
nian is 



^Icad — + + ^LR 



Here 



^L(R) - ^ '^kip)(iipfk{p) 
kip) 

is the Hamiltonian of each of the leads and 

^LR = E^fep^'*4cp + Hc 
kp 



(12) 
(13) 

(14) 



describes the direct coupling between the two leads (pic- 
torially shown as the lower arm of the ring in Fig. [ij . In 
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the phase factor e'* , the flux $ is measured in units of the 
flux quantum, hc/e. Since we use units in which h = 1, 
the flux quantum is lixcj e. Thus, our model Hamiltonian 
is 



U = Hi 



cad 



coup 



(15) 



where the operators cj, cj,, and cj^ (cq, Cj,, and Cp) create 
(annihilate) an electron on the level, on the left lead, and 
on the right lead, respectively, while (6) creates (anni- 
hilates) an excitation of frequency ujq on the molecule. 

In the spirit of the Landaucr approach^ the various 
reservoirs connected to the leads arc described by the 
populations of the excitations pertaining to a specific 
reservoir. The electronic reservoirs on the left and on 
the right of the bridge are characterized by the electronic 
distributions and /j^, 

./l(r)(^) = (l + exp[^L(R)(^ - Ml(r))]) > (16) 

determined by the respective Fermi functions, with 
''^L{R) ~ 1/^B-^L(R)' T^^^ phonon reservoir, which de- 
termines the vibration population on the bridge, is char- 
acterized by the Bose-Einstein distribution. 



with /3y = l/fce^v 



- 1 



(17) 



C. Details of the currents' calculation 

Using the Keldysh technique,— we express the currents 
[see Eqs. ((T5)) and in terms of the Green function 
of the localized level, denoted G = Gq.q(w). The latter 
is calculated up to second order in the coupling to the 
vibrations, 7, 



G< = CMS 



(20) 



Here, (G°) is the retarded (advanced) Green func- 
tion, and is the lesser one (below wc also use the 
greater Green function, G^); Q is the Green function in 
the absence of the coupling with the vibrations. 



(21) 



where EJ" is the self-energy due to the coupling with the 
leads, see Eq. ^ below. Similarly, g< ^ G^'Y^fG''. The 
superscripts on the self energies E have the same meaning 
as those indexing the Green functions. 

We now discuss the two self energies which appear in 
Eq. (PD|) . The first, which comes from the coupling to 
the leads, Ej, is independent of the coupling to the vi- 
brational mode, and is given by 



B. Formal expressions for the currents 

The charge current leaving the left lead is given by 7l = 
—e{NjJ, where N-^ = J^k ^\^k number operator on 

the left reservoir; is the sum of the currents fiowing 
from the left lead to the bridge and into the lower arm 
of the ring, sec Fig. [TJ In terms of the Keldysh Green 
functions!^ 



duj 
2^ 



(18) 



where 



Mu:) = Y,vJg^.^^{u:)-G^.^^{uj) 



5]F,p(e-'*G,,p(c.)-e^*Gp,,(c.) 

kp 



(19) 



The subscripts of the Green functions indicate the corre- 
sponding operators which form them. The electronic en- 
ergy rate in the left electronic reservoir is = {H^}; the 
energy current leaving the left reservoir is given by Eqs. 
(IT5|) and (fTO|) . with e replaced by w (within the integral). 
The electronic currents (charge and energy) associated 
with the right lead are — e(iVp) and —Ej^ = —(Hp); they 
are given by Eqs. ([18]) and (fT9|) upon interchanging k 
with p, L with R, and $ with — <I>. (The explicit fre- 
quency dependence in some of the following equations is 
omitted for brevity.) 



2(1 + A) 



rL + rR-2zVArYr^cos<i> , 



y< 



(/JFl + AFr) + f^{T^ + AFJ 



(1 + A) 

+ 2VAIYrV(/R,-/L)sin<I) 



(22) 



This self energy is determined by the resonance width 
Fl(w) [Fj^(a;)] induced by the coupling with the left 
(right) lead, dressed by all interference paths. In Eqs. 
p2|). A is the (dimcnsionlcss) direct coupling between 
the leads, 

A(c.) = (27r)2^|V^,pP<5(c.-6,)<5(..-ep) . (23) 

k,p 

In our model, the transmission and refiection amplitudes 
of the lower arm of the ring alone (see Fig. [ij, and r^, 
respectively, arc chosen to be real and are given in terms 
of A, 



4A(w) 



(24) 



The total width of the resonance level caused by the 
couplings to the leads, F, is given by — 2ImE[ [sec Eqs. 
^ and (1231)], 



r(c.) = 



1 + A(a;) 



(25) 
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while the asymmetry in those couphngs is characterized 
by 



or alternatively 



[rLM + rj,H]2 



(26) 



(27) 



such that + = 1. None of the above details depends 
on the coupling to the vibrational mode. 

The second self energy in Eq. (^0)) . Sy, comes from 
the coupling with the vibrations. To second order in 7, 
it is given by 

E<(tj) = -f^(^Ng<{uj - Wo) + (1 + N)g<iio + Wo)) 

(28) 



where N denotes the population of the vibrational modes 
[see Eq. (|17|)]: 6ey is the polaron shift, 



J 2n 



(29) 



Note that the Green functions here contain l^f and "Ef" , 
and therefore they also depend on the flux $. 



In terms of these parameters, the integrand determin- 
ing the electronic currents running from left to right [see 
Eqs. (dH) and (gni), and Fig. [T] is 



Jl =(/r - /l) (io[l - rimG°] + t,r,Ta cos $ReG") 

-*^((l + r„a)[G<-/L(G"-G")]+at„sin<i>[G<-/R(G'^-G'-)]) . (30) 

I 



The corresponding integrand determining the currents 
flowing from right to left, J7i?,, is obtained from Eq. ([50)) 
upon interchanging L with R and $ with — a transfor- 
mation which leaves G invariant. 

Let us first consider the sum of the electronic currents, 
flowing from the left towards the ring and from the right 
towards the ring. Using Eqs. (P^ . we find 

+ = (Sr - I]r)G< + S<(G'^ - G^) 

= |G12(S]>S<-S<S>) . (31) 

Therefore, up to second order in the coupling with the 
vibrations, the sum of the currents is given by 

/^(^^ + ^«)^/^i«"i1^?^-^f^?)' 

(32) 

where s = for the sum of the two electronic electric 
currents, and s = 1 for the sum of the two electronic en- 
ergy currents. Inserting here Eqs. ([28]) for the self energy 
due to the coupling with the vibrations, one obtains that 
Eq. p2|l vanishes for s = 0, ensuring the electric current 
conservation. 



On the other hand, for s = 1 one finds 

/ ^(^lM+^rH) -7^-0/ ^|an--)a'^(-+)P 
X (iVS>(c.+ )S<(a._)-(l + 7V)I]>(c._)I]<(c.+)) , 

(34) 

where we have introduced the abbreviations 

a;±=w±^. (35) 

The energy current carried by the vibrations, ujQd{h'^h) jdt 
[see Eq. ([T0| ]. may be calculated from the Keldysli 
phonon Green functions;^ and it is straightforward to 
show that it is given by the result (|34|. ensuring the to- 
tal energy conservation of our model, 

-i^v = K+K- (36) 

D. Explicit expressions for the currents 

Returning to Eq. ([50)1 and inserting there Eqs. (|20p . 
we find that, like the self energies, Jl can be separated 
into two contributions, one from the coupling to the leads 
and the other from the coupling with the vibrations, 



(33) 



Jl(w) ^{^^(yj) ~ h{w)\j\w) + J^{w) , (37) 
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where the 'bare' transition probabihty between the leads 
is (cf. Ref. m 

J\lo) =tl{l- riniG" + ^[l-o? cos2 $]|G°f ) 

+ i„r„racos$ReG"' + ^— IG-'P . (38) 

The first term here, t^, yields the conductance in the ab- 
sence of the ring arm carrying the bridge (see Fig. [T]), 
while the last term yields the conductance of that arm 
alone. All other three terms in Eq. ([55)) result from 
various interference processes which do not involve the 
coupling to the vibrations [apart from their renormaliza- 
tion of the bridge Green function, see Eq. (|20p ]. This 
'bare' transition probability is an even function of the 
flux 

The contribution of the coupling to the vibrations, i.e. 
the last term of Eq. (|37)) , is given to order 7^ by 

J^{^) = -^la-p («+P<(l - /J + E>/J 

+ 5sin$[I]<(l-/i,) + I]>/i,]) , (39) 

where 

a±(a;) = 1 ±S(w)r„(a;) , 
h{u) ^ a{uj)t,{u) , (40) 

see Eqs. ([Ml), and (The function a_ is ob- 

tained from the function a+ upon interchanging L with 

We are now in position to present the detailed expres- 
sions for the currents, by inserting Eq. ([57]) [and the anal- 
ogous one, obtained upon interchanging there L with R 



The various excitation populations characterizing the 
reservoirs are incorporated into the functions F^^, , 

F„„,(^)=7V[l-/„(c^+)]/„,(^_) 

- (l + iV)[l-/„,(c^_ )]/„(..+ ) , (45) 

[see Eqs. ([T5)) and ([T7)) ]. In addition, the four functions 
9±± are given by 

5_l-_l-(a;) = a_j_(tj) ± fe(a;) sin$ (46) 



and 4> with -$] into Eqs. (I6|), d?]), and (HH). The electric 
current is 

-ej'^c{^,^)J^{^), (41) 
while the energy current carried by the electrons is 

^c(c.,$)jV(^). (42) 

The function 

c(..,$)=7^ ^^"+^^^^'"~V ''K)g°(^-)P (43) 

introduced in Eqs. (PT|) and ([^^ has a simple meaning. 
Keeping in mind that (r/2)|t/°p sets the scale for the 
local density of states on the bridge (to leading order), it 
is seen that c(ct', $) is the product of the local densities 
of states at the two shifted frequencies uj± [see Eq. (|35|) ] , 
multiplied by the coupling to the vibrations, 7^. Note 
that the broadening of the energy level on the dot result- 
ing from the coupling to the leads depends on the flux 
$, see Eqs. ([22)) . leading in turn to the flux dependence 
of Q and the function c, Eq. P5|) . This is the cos$ de- 
pendence resulting from the usual interference processes, 
discussed at the end of Sec. [ill 

The other function introduced above is , where s — 
in Eq. (gT]), and s = 1 in Eq. (gH), 



(44) 



for the equal-sign functions g, and 

.g±^(w) = a_,_(cL;) =F ^(w) sin$ (47) 

for the unequal ones. 

Finally, the energy current carried by the vibrations is 



+ ^RL (^) ( -^+.9— (^+ ) - 1'^- ).9-+ (w+ ) 
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given by 



in 



+ ^'RR(w).g„+(w+)ff_+(w_) 
+ ^lR(w).g+„(w+)5_+(a;_) 



(48) 



These lengthy expressions are significantly simplified 
when one assumes a symmetrically-coupled bridge, or a 
perfect bare transmission of the lower arm of the ring 
(see Fig. [1]). In both cases, Eqs. pOj) yield a± — 1. 



IV. THE LINEAR-RESPONSE REGIME 

In the linear-response regime one expands all cur- 
rents to first order in the driving forces, which in our 
case are A/x, AT, and ATy (see discussion in Sec. 
HI]) . This amounts to expanding the Fermi distribu- 
tions, Eqs. (|16p . around the common temperature T 
and the common chemical potential ^ of the device, i.e., 
around /(w) = [exp(/3(aj — and the vibrational 
mode population, Eq. (ITTl) . around the Bose function 
= [exp(/3wo) - where ^ = 

The electric current in the linear-response regime takes 
the form 



7 = G($)^+K($)^+xV(ci>)^ 
el 1 



(49) 



As is clear from the discussion in Sec. IIIIl there arc two 
contributions to each of the first two terms of the electric 
current, one resulting from the coupling to the leads [the 
first term on the right-hand side of Eq. (HI])] and the 
other coming from the coupling to the vibrations. The 
last term in Eq. (|49p arises from the latter coupling alone. 
Accordingly, 



with 



where 



G = G' + , 
K = K' -(- , 



(50) 



(51) 



(52) 



is essentially the transition probability, Eq. (|38)) . at the 
Fermi energy (at low temperatures). Similarly the trans- 
port coefficient related to the thcrmopower is 



K'($) 



2^ 



(53) 



These expressions for the transport coefficients are the 
standard ones, see for example Ref. [Ill; both coefficients, 
G' and K', are even in the magnetic flux. 



Next we examine the contribution of the 'vibronic' 
transitions to the transport coefficients. Introducing the 
function 

r^(c^, $) = /3iVT/(w_)[l - /(w+)]c(w, $) , (54) 
we find 

G^($) = j ^r^iuj, $)mo(c^, $) , (55) 

where 

+ a(w_)a(w+)t„(a;_)t„(cj_^)sin2$ . (56) 

[We remind the reader that a and a characterize the 
asymmetry of the bridge, see Eqs. ([?7)) . while and 
arc the transmission and reflection amplitudes of the 
lower arm of the ring, see Eqs. (j24p . The frequencies ui± 
are given in Eq. ([35]).] Notice that the full conductance 
G'-|-G^ is even in the flux. Thus it is seen that the 
electron-vibration interaction introduces only a modest 
modification in the conductance, but does not lead to 
major effects. 

This is not the case with all other transport coefficients 
of the ring. The thcrmopower coefficient is given by 

KV(ci,) = e J ^(c.-/i)r^(w,$)mo(w,$) 

+ e^sin$ / — r^(w,$)mi(cj) , (57) 
2 J -K 

where 

-a(w„)to(w_)a(a;+)r„(w+) . (58) 

While the first term of is even in the fiux, the second 
term is an odd function of $. In order not to vanish, this 
odd component necessitates that the bridge will not be 
coupled symmetrically to the leads, and that the trans- 
mission of the lower arm of the ring will not be perfect 
(to ensure interference). In addition, is invariant un- 
der the transformation LoR [which changes the sign of 
a, see Eq. dST])] and $ 4-> 

Finally, the coefficient [see Eq. (j49])] is 



r'{<S>) = eoj,J ^r^(c.,$)TO2(c.,<i>) , 



(59) 



with 



TO2(<^7 ^) ~ Q;('^-)''o('^-) + Q^('^-)^o('^-) sin$ 

- a(w+)r<,(w+) - a(a;+)t„(cj+)sin$ . (60) 

Note that under the transformation LoR and $<->—<!>, 
changes its sign. Also note that, again, the sin$ 
dependence of this transport coefficient necessitates the 
coupling to the vibrational modes. 
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In an analogous fashion we derive the heat current car- 
ried by the electrons, and find 

/q-K(-$)^+K2($)— +X^($)^ . (61) 

Here, the transport coefficient related to the thermal con- 
ductivity, K2, has again two contributions. 



EXAMPLES AND DISCUSSION 



K2 = + Ky , 



(62) 



where 



K^(*)= / ^(c.-m)V(c.,$) (63) 
as expected, while the 'vibronic' term is 

Ky($)^ / ^(c.-M)'r^(^,$)mo(c.,$) 

+ t/ v'^^(-.'i>)-o(-,*), (64) 

with 

- a{uj_)a{u;M^M^+) sin' * ■ (65) 

One notes that K2 is even in the flux. The last transport 
coefficient in Eq. (|6T|) is 



T^{u,<i>)m^{u,^) , (66) 



with 



m2(a;, $) = a{u!_)rg{u!_) + a{u!_)t^{uj_) sin $ 

+ a(w+)r„(^+)+a(w+)t„(^+)sin$ . (67) 

This coefficient changes its sign under the transforma- 
tion L4->R and $ — $, as expected from the Onsagcr 
relations. 

The energy current from the thermal bath takes the 
form 



i?v=xV(-*)— +X^(-$ 



AT „w AT, 



C 



e T T 

where the third coefficient here is 



, (68) 



;/^r>,<i>) 



CV^2c.2 / (,.,$) , (69) 



and is an even function of the flux. Collecting all results, 
we arrive at the matrix form for the relations among the 
currents and the driving forces, Eqs. ([5]) and 



As is clear from the results in Sec. IIVI (and also from 
Rcf. [H, which treated the bridge without the lower 
branch), the interesting effect induced exclusively by the 
coupling of the electrons to the vibrational modes is the 
possibility to create an electric current, or an electronic 
heat current, by applying a temperature difference ATy 
on the phonon bath thermalizing this mode. These new 
thermoelectric phenomena are specifled by the two co- 
efficients and X^. Eqs. ([59|) and (|66p . respectively. 
All other transport coefficients related to the electronic 
currents are mainly due to the transport of the electrons 
between the electronic terminals, with slight modifica- 
tions from the (small) coupling to the vibrations. We 
therefore confine the main discussion in this section to 
the coefficients X^ and X^. To make a closer connection 
with a possible experiment, we introduce the (dimension- 
less) coefficients 



= e/3- 



.XV 
G 



and 



?v - — 



(70) 



(71) 



The first gives the potential drop across the molecular 
bridge created by ATy when the temperature drop there, 
AT, vanishes, and the second yields the temperature dif- 
ference created by ATy when A^ = [or the inverse pro- 
cesses, see Eqs. ([5]) and (O]. For both the conductance, 
G, and the thermal conductance, K2, we use below their 
leading terms, resulting from the coupling to the leads 
alone (the numerators result from the coupling to the 
vibrations, and hence are already of order 7^). 

As mentioned above, the transport coefficients of our 
three-terminal junction obey the Onsager-Casimir rela- 
tions. They do it though in a somewhat unique way: 
the "off-diagonal" elements arc related to one another by 
the reversal of the magnetic field. However, they are not 
a purely odd function of it. A special situation arises 
when the molecule is connected symmetrically to the two 
leads. In that case, the anisotropy parameter a vanishes, 
while a = 1 [sec Eqs. ^ and Then, the two 

coefficients, X^ and X^, are odd functions of the flux, 
resembling the thermal Hall effect discussed recently in 
connection with quantum magnets,— 

XV(<i>)=ec.oSin$ J ^T"" {u,^){t,{uj^) - t,{iu^)) , 

(72) 



and 



X^($)=WoSin$ / — r^(w,$) 



{t,{uj_)iuj_ - ^) - <„(^+)(w+ - fi)) . (73) 
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In other words, the thermoelectric processes described by 
X^($) and X^($) require a certain symmetry-breaking. 
In the absence of the magnetic field, that is supplied by 
the spatial asymmetry of the junction; in the presence of 
a flux, those processes appear also for a junction symmet- 
rically coupled to the leads, provided that the couplings 
to the leads depend on the energy. 

When the two leads connected to the electronic reser- 
voirs are identical (making the molecular bridge sym- 
metric) the transmission amplitude of the direct bond 
between the two leads, [see Eq. ([24)) ]. is an even func- 
tion of w. The function [uj, $) given by Eqs. P5|) and 
(|54p . is not entirely even or odd in uj, and therefore a 
priori the integrals which give X^($) and X^($) do not 
vanish. However, the asymmetry in the w— dependence 
of the integrand (which results from the w— dependence 
of the Green function) is not significant. As a result, 
is extremely small, while is not (because of the 
extra w factor in the integrand), see Fig. [2] These 
plots are computed using r(a;) = rVl-(w/M^)2, and 
A(w) = A"[l - {io/W)\ where W is half the bandwidth, 
and all energies are measured in units of (3 = l/(fcgr) 
(we have set T" = A'' = 1 and = 50). 



5.x 10"'' 




FIG. 2: (color online.) The transport coefficients and 
as functions of the flux (measured in units of the flux 
quantum) and Puiq, for a symmetric bridge. 

The relative magnitudes of and are significantly 
changed when the molecule is coupled asymmetrically to 
the leads, and moreover, the character of the charge car- 



riers on the two reservoirs is different. Let us assume 
that the left reservoir is represented by an electron band, 
such that the partial width it causes to the resonant level 
is given by 



(74) 



while the right reservoir is modeled by a hole band, with 



rR,(^) = ro. 



(75) 



The corresponding quantity pertaining to the lower arm 
of the ring in Fig. [T] is 





— UJ 


1. 








1 





(76) 



Here, uj^ is the bottom of the conductance band (on the 
left side of the junction), while is the top of the hole 
band (on the right one). The energy integration deter- 
mining the various transport coefficients is therefore lim- 
ited to the region w^, < w < w„. (For convenience, we 
normalize the F's by the full band width, a;„ — uj^.) Ex- 
emplifying results in such a case are shown in Fig. [31 
computed with F^ = = A° = 1 [see Eqs. (d, ([75]), 
and (j76p ]. and cj^ = — w,„ = 100, all in units of (3. 

S^KPyf 




FIG. 3: (color online.) Same as Fig. [2l for an asymmetric 
bridge. 

Remarkably enough, the coefficient S^, which mea- 
sures the ability of the device to turn a temperature 
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difference on the phonon bath thermalizing the molec- 
ular vibrations into a temperature difference across the 
molecule is not very sensitive to the details of the model, 
apparently because the term proportional to sin <J> in Eq. 
(|66p is the dominant one. This means that by reversing 
the direction of the magnetic field, one reverses the sign 
of the temperature difference or alternatively, the sign of 
the electronic heat current. On the other hand, the coef- 
ficient S^, which sets the scale of the capability to turn a 
temperature difference across that phonon bath into an 
electric current, is far more sensitive to the details of the 
model (as expressed e.g. in our choice for the density of 
states on the leads) and is less affected by the magnetic 
field. 

We now return to the diagonal charge and heat con- 
ductances. Figure HI shows the flux-dependence of their 
dimcnsionlcss ratio, /3^e^K2/G, plotted with the same 
parameters as above. It is interesting to observe that, for 
both the symmetric and the asymmetric bridges, this ra- 
tio remains quite close to the textbook Wiedemann- Franz 
ratio 7r^/3, calculated from the Fermi-Dirac distribution 
of a free electron gas4^ This is quite surprising, in view of 



the much richer resonance structure of the transmission 
through the ring. 

In conclusion, we have found that the themoelectric 
transport coefficients through a vibrating molecular junc- 
tion, placed on an Aharonov-Bohm interferometer, have 
an interesting dependence on the magnetic fiux. In par- 
ticular, the coefficients which relate the temperature dif- 
ference between the phonon and electron reservoirs to the 
charge and heat currents carried by the electrons, which 
exist only due to the electron-vibron coupling, can be 
enhanced by the magnetic flux. 
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